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Abstract 

Given a hypergraph G, we introduce a Grassmann algebra over the ver- 
tex set, and show that a class of Grassmann integrals permits an expansion in 
terms of spanning hyperforests. Special cases provide the generating functions 
for rooted and unrooted spanning (hyper)forests and spanning (hyper)trees. 
All these results are generalizations of Kirchhoff 's matrix-tree theorem. Fur- 
thermore, we show that the class of integrals describing unrooted spanning 
(hyper)forests is induced by a theory with an underlying 0SP(1|2) supersym- 
metry. 
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Keywords: Graph, hypergraph, forest, hyperforest, matrix-tree theorem, Grass- 
mann algebra, Grassmann integral, cr-model, supersymmetry. 



1 



1 Introduction 



Kirchhoff's matrix-tree theorem [1-3] and its generalizations [4-6], which express 
the generating polynomials of spanning trees and rooted spanning forests in a graph 
as determinants associated to the graph's Laplacian matrix, play a central role in 
electrical circuit theory [7, 8] and in certain exactly-soluble models in statistical me- 
chanics [9, 10]. 

Like all determinants, those arising in Kirchhoff's theorem can be rewritten as 
Gaussian integrals over fermionic (Grassmann) variables. Indeed, the use of Grassmann- 
Berezin calculus [11] has provided an interesting short-cut toward the classical matrix- 
tree result as well as generalizations thereof [6,12]. For instance, Abdesselam [6] has 
obtained in a simple way the recent pfaffian-tree theorem [13-15] and has generalized 
it to a hyperpfaffian-cactus theorem. 

In a recent letter [12] we proved a far-reaching generalization of Kirchhoff's the- 
orem, in which a large class of combinatorial objects are represented by suitable 
non-Gaussian Grassmann integrals. In particular, we showed how the generating 
function of spanning forests in a graph, which arises as the g — > limit of the partition 
function of the g-state Potts model [16-19], can be represented as a Grassmann inte- 
gral involving a quadratic (Gaussian) term together with a special nearest-neighbor 
four-fermion interaction. Furthermore, this fermionic model possesses an 0SP(1|2) 
supersymmetry. 

This fermionic formulation is also well-suited to the use of standard field-theoretic 
machinery. For example, in [12] we obtained the renormalizat ion-group flow near 
the spanning-tree (free-field) fixed point for the spanning-forest model on the square 
lattice, and in [20] this was extended to the triangular lattice. 

In the present paper we would like to extend the fermionic representation of 
spanning forests from graphs to hypergraphs. Hypergraphs are a generalization of 
graphs in which the edges (now called hyperedges) can connect more than two ver- 
tices [21-23]. In physics, hypergraphs arise quite naturally whenever one studies a 
fc-body interaction with > 2|j We shall show here how the generating function of 
spanning hyperforests in a hypergraph, which arises as the g — > limit of the partition 
function of the g-state Potts model on the hypergraph [24], can be represented as a 
Grassmann integral involving a quadratic term together with special multi-fermion 
interactions associated to the hyperedges. Once again, this fermionic model possesses 
an 0SP(1|2) supersymmetry. This extension from graphs to hypergraphs is thus not 
only natural, but actually sheds light on the underlying supersymmetry. 

Let us begin by recalling briefly the combinatorial identities proven in [12], which 
come in several levels of generality. Let G = (V, E) be a finite undirected graph with 
vertex set V and edge set E. To each edge e we associate a weight We, which can 



^ For examples in the recent physics literature where the hypergraph concept is used, see for 
instance [24-27]. 
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be a real or complex number or, more generally, a formal algebraic variable; we then 
define the (weighted) Laplacian matrix L = (Ljj)jjgy for the graph G by 



(1.1) 



We introduce, at each vertex i E V, a pair of Grassmann variables ipi, ipi, which 
obey the usual rules for Grassmann integration [11,28]. Our identities show that 
certain Grassmann integrals over and ip can be interpreted as generating functions 
for certain classes of combinatorial objects on G. 

Our most general identity concerns the operators Qr associated to arbitrary con- 
nected subgraphs T = (yr,Er) of G via the formula 



Qr 



n 



n 



:i.2) 



(Note that each Qr is even and hence commutes with the entire Grassmann algebra.) 
We prove the very general identity 



exp 



ijLijj + ^trQr 



J? spanning CG VeG-ff / a=l 



where the sum runs over spanning subgraphs H C G consisting of connected compo- 
nents {Hi, . . . , Hi), and the weights W{Ha) are defined by 



W{H, 



:i.4) 



where F -< Ha means that Ha contains F and contains no cycles other than those 
lying entirely within F. 

Let us now specialize (11. 3p to the case in which tr = ti when F consists of a single 
vertex i with no edges, tr = Ue when F consists of a pair of vertices i,j linked by an 
edge e, and = otherwise. We then have 



T>{ip,ip) exp ipL^p + "^titpitpi + UijWij^i^jipjipj 



E n »« n E E 



Uf, 



[1.5) 



e£E{Fa) 
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where the sum runs over spanning forests F in G with components Fi, . . . , F^; here 
V{Fa) and E{Fa) are, respectively, the vertex and edge sets of the tree Fa. 

If we further speciahze (11.51) to Ue = —A for all edges e (where A is a global 
parameter), we obtain 

V{ip,ip) exp 4}Li; + Utpj^j - X^Wiji^i^ii^j^j 



;i.6) 



F G :f{G) 

F = {Fi,...,Fi) 



z£F 



a=l 



ieV(F^) 



since \E{Fa) \ = \V{Fa) \ — 1. If, in addition, we take = A for all vertices i, then we 
obtain 

'D{ijj,tp) exp ^pL^p + X'^tpitpi - A Wij^jipjipjipj 



= E 

F&r{G) \e&F 

= ^'^'E (n 

F&T{G) \e&F 



XHF) 



We 

X 



[1.7a) 



;i.7b) 



where k{F) is the number of component trees in the forest F; this is the generating 
function of (unrooted) spanning forests of G. Furthermore, as discussed in [12] and 
in more detail in Section [7] below, the model (ll.7p possesses an 0SP(1|2) invariance. 
If, by contrast, in (II. 6p we take A = but {tj} general, we obtain 



V{il),%l)) exp ipLip + tjipi^i 



E n-. n E 



(1.8) 



F 6 T{G) 
F = {Fi,...,Ft) 



26F 



which is the formula representing rooted spanning forests (with a weight ti for each 
root i) as a fermionic Gaussian integral (i.e., a determinant) involving the Lapla- 
cian matrix (this formula is a variant of the so-called "principal-minors matrix-tree 
theorem" ) . 

In this paper we shall not attempt to find the hypergraph analogue of the general 
formula (11.31) . but shall limit ourselves to finding analogues of (ll.5p -( |T78l) . The formu- 
lae to be presented here thus express the generating functions of unrooted or rooted 
spanning hyperforests in a hypergraph in terms of Grassmann integrals. In particular, 
the hypergraph generalization of (ll.7p possesses the same 0SP(1|2) supersymmetry 
that (I1.7P does (see Section [T]). 
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The proof given here of all these identities is purely algebraic (and astonishingly 
simple); the crucial ingredient is to recognize the role and the rules of a certain 
Grassmann subalgebra (see Section HI). It turns out (Section [7]) that this subalgebra 
is nothing other than the algebra of 0SP(l|2)-invariant functions, though this is far 
from obvious at first sight. The unusual properties of this subalgebra (see Lemma SH]) 
thus provide a deeper insight into the identities derived in [12] as well as their gen- 
eralizations to hypergraphs, and indeed provide an alternate proof of f ll.6l) -( fL8l) . 
Pictorially, we can say that it is the underlying supersymmetry that is responsible for 
the cancellation of the cycles in the generating function, leaving only those spanning 
(hyper)graphs that have no cycles, namely, the (hyper) forests. 

In particular, the limit of spanning hyper irees, which is easily extracted from the 
general expression for (rooted or unrooted) hyperforests, corresponds in the 0SP(1|2)- 
invariant cx-model to the limit in which the radius of the supersphere tends to infinity, 
so that the nonlinearity due to the curvature of the supersphere disappears. However, 
the action is in general still non-quadratic, so that the model is not exactly soluble. 
(This is no accident: even the problem of determining whether there exists a spanning 
hypertree in a given hypergraph is NP-complete [29].) Only in the special case of 
ordinary graphs is the action purely quadratic, so that the partition function is given 
by a determinant, corresponding to the statement of Kirchhoff 's matrix-tree theorem. 

The 0SP(l|2)-invariant fermionic models discussed in [12] and the present paper 
can be written in three equivalent ways: 

• As purely fermionic models, in which the supersymmetry is somewhat hidden. 

• As (T-models with spins taking values in the unit supersphere in ]Ri|2, in which 
the supersymmetry is manifest. 

• As A^- vector models [= 0(A^)-symmetric cr-models with spins taking values in 
the unit sphere of M^] analytically continued to A^ = — 1. 

The first two formulations (and their equivalence) are discussed in Section [7] of the 
present paper. Further aspects of this equivalence — notably, the role played by the 
Ising variables arising in (17.91) and neglected here — will be discussed in more detail 
elsewhere [30]. 

In a subsequent paper [31] we will discuss the Ward identities associated to the 
0SP(1|2) supersymmetry, and their relation to the combinatorial identities describing 
the possible connection patterns among the (hyper)trees of a (hyper)forest. 

The method proposed in the present paper has additional applications not con- 
sidered here. With a small further effort, a class of Grassmann integrals wider than 
f l5.2p / (16.31) — allowing products Yl^ f^^ in the action in place of the single operators 

— can be handled. Once again one obtains a graphical expansion in terms of 
spanning hyperforests, where now the weights have a more complicated dependence 
on the set of hyperedges, thus permitting a description of certain natural interaction 
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patterns among the hyperedges of a hyperforest (see Remark at the end of Section [3]). 
This extended model is, in fact, the most general Hamihonian that is invariant under 
the 0SP(1|2) supersymmetry. 

The plan of this paper is as follows: In Section [5] we recall the basic facts about 
graphs and hypergraphs that will be needed in the sequel. In Section [3] we define the 
g-state Potts model on a hypergraph and prove the corresponding Fortuin-Kasteleyn 
representation. (This section is unnecessary for the proof of the combinatorial iden- 
tities that form the main focus of this paper, but it provides additional physical 
motivation.) In Section H] we introduce the Grassmann algebra over the vertex set V, 
and study a subalgebra with interesting and unusual properties, which is generated 
by a particular family of even elements with A ^ V. In Section [5] we study a 
very general partition function involving the operators and we show how it can 
be expressed as a generating function of spanning hyperforests in a hypergraph with 
vertex set V. In Section [6] we study a somewhat more general Grassmann integral, 
which can be interpreted as a correlation function in this same Grassmann model; we 
show how it too can be expressed as a sum over spanning hyperforests. In Section [7] 
we show that in one special case — namely, the hypergraph generalization of (11.71) — 
the model studied in the preceding sections can be rewritten as an 0SP(l|2)-invariant 
cr-model, and indeed is the most general 0SP(l|2)-invariant Hamiltonian involving 
ip and ip. These cr-model formulae motivate the definition of f^^^ given in Section IU 
which might otherwise remain totally mysterious. 

In Appendix 1X1 we prove a determinantal formula for In Appendix [B] we 

present a graphical formalism for proving both the classical matrix-tree theorem and 
numerous extensions thereof, which can serve as an alternative to the algebraic ap- 
proach used in the main body of this paper. 

Let us stress that everything in this paper is mathematically rigorous, with the 
possible exception of Section [71 Mathematicians unfamiliar with the Grassmann- 
Berezin calculus can find a brief introduction in [6, Section 2] or [32, Appendix A]. 

2 Graphs and hypergraphs 

A (simple undirected finite) graph is a pair G = {V, E) , where y is a finite set 
and E IS a. collection (possibly empty) of 2-element subsets of V^H The elements of 
V are the vertices of the graph G, and the elements of E are the edges. Usually, in 
a picture of a graph, vertices are drawn as dots and edges as lines (or arcs). Please 
note that, in the present definition, loops ( <Z) ) and multiple edges ( <II>) are not 

^ To avoid notational ambiguities it should also be assumed that E DV = 0. This stipulation is 
needed as protection against the mad set theorist who, when asked to produce a graph with vertex 
set V = {0, 1, 2}, interprets this a la von Neumann as = { 0, {0}, {0, {0}} }, so that the vertex 2 
is indistinguishable from the edge {0, 1}. 
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allowedH We write \V\ (resp. \E\) for the cardinahty of the vertex (resp. edge) set; 
more generally, we write 15*1 for the cardinality of any finite set 5*. 

A graph G' = {V, E') is said to be a subgraph of G (written G' C G) in case 
V' '^V and E' C E. If V = V, the subgraph is said to be spanning. We can, by a 
slight abuse of language, identify a spanning subgraph (V, E') with its edge set E'. 

A walk (of length k > 0) connecting vq with Vk in G is a sequence {vq, ei, f i, 62, ^2, 
. . . , ek,Vk) such that all Vi G V, all Cj G E, and Vi-i,Vi E Ci for 1 < i < k. A path 
in G is a walk in which vq, . . . ,Vk are distinct vertices of G and ei, . . . , Cfc are distinct 
edges of G. A cycle in G is a walk in which 

(a) f 0, . . . , f fe_i are distinct vertices of G, and Vk = f 

(b) ei, . . . , Cfc are distinct edges of G; and 

(c) > 20 

The graph G is said to be connected if every pair of vertices in G can be connected 
by a walk. The connected components of G are the maximal connected subgraphs of G. 
It is not hard to see that the property of being connected by a walk is an equivalence 
relation on V, and that the equivalence classes for this relation are nothing other 
than the vertex sets of the connected components of G. Furthermore, the connected 
components of G are the induced subgraphs of G on these vertex setsE] We denote 
by k{G) the number of connected components of G. Thus, k{G) = 1 if and only if G 
is connected. 

A forest is a graph that contains no cycles. A tree is a connected forest. (Thus, 
the connected components of a forest are trees.) It is easy to prove, by induction on 
the number of edges, that 

\E\ - \V\ + k{G) > (2.1) 

for all graphs, with equality if and only if G is a forest. 

In a graph G, a spanning forest (resp. spanning tree) is simply a spanning subgraph 
that is a forest (resp. a tree). We denote by J-'{G) [resp. T{G)] the set of spanning 
forests (resp. spanning trees) in G. As mentioned earlier, we will frequently identify 
a spanning forest or tree with its edge set. 



^ This restriction is made mainly for notational simplicity. It would be easy conceptually to allow 
multiple edges, by defining E as a multiset (rather than a set) of 2-element subsets of V (cf. also 
footnote [7] below) . 

Actually, in a graph as we have defined it, all cycles have length > 3 (because ei ^ 62 and 
multiple edges are not allowed). We have presented the definition in this way with an eye to the 
corresponding definition for hypergraphs (see below), in which cycles of length 2 are possible. 

^ If T^' C V, the induced subgraph of G on V, denoted G'[y], is defined to be the graph {V , E') 
where E' is the set of all the edges e G i? that satisfy e C y (i.e., whose endpoints are both in V'). 
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Figure 1: A forest (left) and a hyperforest (right), each with four components. Hy- 
peredges with more than two vertices are represented pictorially as star-hke polygons. 



Finally, we call a graph unicyclic if it contains precisely one cycle (modulo cyclic 
permutations and inversions of the sequence vo,ei, vi,e2,V2, ■■■ ,ek,Vk)- It is easily 
seen that a connected unicyclic graph consists of a single cycle together with trees 
(possibly reduced to a single vertex) rooted at the vertices of the cycle. 

Hypergraphs are the generalization of graphs in which edges are allowed to con- 
tain more than two vertices. Unfortunately, the terminology for hypergraphs varies 
substantially from author to author, so it is important to be precise about our own 
usage. For us, a hypergraph is a pair G = (V, E), where is a finite set and E is a 
collection (possibly empty) of subsets of V, each of cardinality >2^ The elements of 
V are the vertices of the hypergraph G, and the elements of E are the hyperedges (the 
prefix "hyper" can be omitted for brevity). Note that we forbid hyperedges of or 1 
vertices (some other authors allow these) III We shall say that A E E is a k-hyperedge 
if y4 is a /c-element subset of V. A hypergraph is called k-uniform if all its hyperedges 
are /c-hyper edges. Thus, a graph is nothing other than a 2-uniform hypergraph. 

The definitions of subgraphs, walks, cycles, connected components, trees, forests 
and unicyclics given above for graphs were explicitly chosen in order to immediately 
generalize to hypergraphs: it suffices to copy the definitions verbatim, inserting the 
prefix "hyper" as necessary. See Figure [1] for examples of a forest and a hyperforest. 

The analogue of the inequality fl2.ll) is the following: 



^ To avoid notational ambiguities it is assumed once again that E OV = 0. 

^ Our definition of hypergraph is the same as that of McCammond and Meier [33] . It is also the 
same as that of Grimmett [24] and Gessel and KaHkow [34] , except that these authors allow multiple 
edges and we do not: for them, E is a multiset of subsets of V (allowing repetitions), while for us 
E is a set of subsets of V (forbidding repetitions). 
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Proposition 2.1 Let G = {V,E) be a hypergraph. Then 

- \V\ + k{G) > 0, (2.2) 

with equality if and only if G is a hyperforest. 

Proofs can be found, for instance, in [22, p. 392, Proposition 4] or [34, pp. 278-279, 
Lemma] . 

Please note one important difference between graphs and liypergraphs: every 
connected graph has a spanning tree, but not every connected hypergraph has a 
spanning hypertree. Indeed, it follows from Proposition 12.11 that if G is a A;-uniform 
connected hypergraph with n vertices, then G can have a spanning hypertree only if 
k — 1 divides n — 1. Of course, this is merely a necessary condition, not a sufficient 
one! In fact, the problem of determining whether there exists a spanning hypertree 
in a given connected hypergraph is NP-complete (hence computationally difficult), 
even when restricted to the following two classes of hypergraphs: 

(a) hypergraphs that are linear (each pair of edges intersect in at most one vertex) 
and regular of degree 3 (each vertex belongs to exactly three hyperedges); or 

(b) 4-uniform hypergraphs containing a vertex which belongs to all hyperedges, and 
in which all other vertices have degree at most 3 (i.e., belong to at most three 
hyperedges) 

(see [29, Theorems 3 and 4]). It seems to be an open question whether the problem 
remains NP-complete for 3-uniform hypergraphs. 

Finally, let us discuss how a connected hypergraph can be built up one edge at a 
time. Observe first that if G = (V, E) is a hypergraph without isolated vertices, then 
every vertex belongs to at least one edge (that is what "without isolated vertices" 

means!), so that V = [J A. In particular this holds if G is a connected hypergraph 

AeE 

with at least two vertices. So let G = (V, E) be a connected hypergraph with \ V\ > 2; 

let us then say that an ordering {Ai, . . . , Am) of the hyperedge set E is a. construction 

e 

sequence in case all of the hypergraphs G^ = ( |J Ai, {Ai, . . . , Ae}) are connected 

1=1 

(1 < £ < m). An equivalent condition is that ( [J Aij H Ai ^ for 2 < i < m. We 

^ i=i ^ 

then have the following easy result: 

Proposition 2.2 Let G = {V, E) he a connected hypergraph with at least two vertices. 
Then: 

(a) There exists at least one construction sequence. 
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(b) If G is a hypertree, then for any construction sequence [Ai,...,Am) we have 

[ [J AAnAe =1 for all£ {2<e <m). 
^ i=i ' 

(c) IfG is not a hypertree, then for any construction sequence {Ai, . . . , Am) we have 



( [JAAnA, 

^ i=i ' 



> 2 for at least one i. 



Proof, (a) The "greedy algorithm" works: Let Ai be any hyperedge; and at each 
stage i > 2, let Ai be any hyperedge satisfying I [J Ai) CiAi ^ (such a hyperedge 



i=l 



has to exist, or else G fails to be connected). 

(b) and (c) are then easy consequences of Proposition 12.11 □ 



3 Potts model on a hypergraph 

Let g be a positive integer, and let S' be a set of cardinality q. Then the q-state 
Potts model on the hypergraph G = {V,E) is defined as follows [24]: At each vertex 
i & V we place a color (or spin) variable cr.j G S. These variables interact via the 
Hamiltonian 

^Potts(cr) = - ^ JA^Aif^) , (3.1) 

where {Ja}a£E are a set of couplings associated to the hyperedges of G, and the 
Kronecker delta 6a is defined for A = {ii, . . . , i^} by 

= I' (3.2) 

I U otherwise 

The partition function Z^"***^ is then the sum of exp[— 7ipotts(o")] over all configurations 
o" = {ai}i(zv. 

It is convenient to introduce the quantities va = e"^^ — 1; we write v = {va}a<^e 
for the collection of hyperedge weights. We can then prove the Fortuin-Kasteleyn 
(FK) representation [35,36] for the hypergraph Potts model [24], by following exactly 
the same method as is used for graphs (see e.g. [19, Section 2.2]): 

Proposition 3.1 (Fortuin-Kasteleyn representation) Let G = {V, E) be a hy- 
pergraph. Then, for integer q > 1, we have 



cr: V^S E'CE A&E' 

where k{E') denotes the number of connected components in the hypergraph {V,E'). 
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Proof. We start by writing 

^S°"^(9,v) = 5^ exp[-Hpotts(a)] = Yl n + ■ (3.4) 

Now expand out the product over A & E, and let E' (1 E he the set of hyperedges for 
which the term VaSa{o') is taken. Now perform the sum over configurations {crjjjgy: 
in each connected component of the spanning subhypergraph (V, E') the color aj must 
be constant, and there are no other constraints. Therefore, 

^r^(g,v) = E^'^""'^ n^^' (3-5) 

E'CE A£E' 

as was to be proved. □ 

Please note that the right-hand side of (13.31) is a polynomial in g; in particular, 
we can take it as the definition of the Potts- model partition function ZG{q,\) for 
noninteger q. 

Let us discuss in particular the various types of g ^ limits that can be taken 
in the hypergraph Potts model, by following a straightforward generalization of the 
method that is used for graphs [19, Section 2.3]: 

The simplest limit is to take q ^ with fixed v. From the definition fl3.3p we 
see that this selects out the spanning subhyper graphs E' ^ E having the smallest 
possible number of connected components; the minimum achievable value is of course 
k{G) itself (= 1 in case G is connected, as it usually is). We therefore have 

limg"'^(^)ZG(g,v) = Cg(v) , (3.6) 

where 

^g(v) = e n (3-7) 

E' CE A&E' 
k{E') = k{G) 

is the generating polynomial of "maximally connected spanning subhypergraphs" (= 
connected spanning subhypergraphs in case G is connected). 

A different limit can be obtained by taking g — with fixed values of wa = 
VA/q^^^~^- From (13. 3p we have 

k{E')+ J2 (1^1-1) -r-r 

^G(g,{gi^i-'ti^A}) = 5Z g JJwa (3.8) 

E'CE AeE' 

Using now Proposition 12.11 we see that the limit g — selects out the spanning 
hyperforests: 

limg-I^IZG(g,{gl^l-^«;A}) = Fg(w) , (3.9) 
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where 

E n^^ (3.10) 

E'£J^(G) A&E' 

is the generating polynomial of spanning hyperforests . 

By a further limit we can obtain spanning hypertrees. To see this, assume first that 
G is connected (otherwise there are no spanning hypertrees). In Cg(v), replace va 
by X''^^~^va and let A — > 0; then we pick out the connected spanning subhypergraphs 
having the minimum value of X]a6-B'(I^I ~^)^ which by Proposition 12. II are precisely 
the spanning hypertrees: 

limA'(l^l-i)CG({Al^l-V}) = Tci^r), (3.11) 

A— >0 

where 

TgM = E n (3.12) 

E'€T{G) A£E' 

is the generating polynomial of spanning hypertrees . Alternatively, in Fg(w), replace 
Wa by X^^^^^wa and let A oo; then we pick out the spanning hyperforests having 
the maximum value of X^^g£;/(|^| — 1), which by Proposition 12.11 are those with the 
minimum number of connected components, i.e. again the spanning hypertrees: 

hm \-^\''\-^^Fa{{X\^\~'wA]) = Tg{w) . (3.13) 

A— >oo 

There is, however, one important difference between the graph case and the hyper- 
graph discussed in Section [21 every connected graph has a spanning tree, but 
not every connected hypergraph has a spanning hypertree. So the limits (13. lip and 
f l3.13p can be zero. 



4 A Grassmann subalgebra with unusual proper- 
ties 

Let be a finite set of cardinality n. For each i E V we introduce a pair ipi, ipi of 
generators of a Grassmann algebra (with coefficients in R or C). We therefore have 
2n generators, and the Grassmann algebra (considered as a vector space over R or 
C) is of dimension 2'^^. 

For each subset A C ^, we associate the monomial ta = YlieA i'ii^ij where = 1. 
Please note that all these monomials are even elements of the Grassmann algebra; in 
particular, they commute with the whole Grassmann algebra. Clearly, the elements 
{'Ta}acv span a vector space of dimension 2"'. In fact, this vector space is a subalgebra, 
by virtue of the obvious relations 

(taub if Ar]B = 

TATE = { (4.1) 

lO ifAnfi^0 
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Let us now introduce another family of even elements of the Grassmann algebra, 
also indexed by subsets of V, which possesses very interesting and unusual properties. 
For each subset A ^ V and each number A (in M or C), we define the Grassmann 
element 

/i") = A(l-|A|)r^ + ^r^,, - J2 M-^Axfe-}- (4.2) 

i + j 

(The motivation for this curious formula will be explained in Section [71) For instance, 
we have 

= A (4.3a) 
ffl = 1 foralH (4.3b) 

= -XiJii)iijji)j + ijj){ipi- ipj) (4.3c) 

and in general 



"=1 l<a,/9<fc 

(4.4) 

(Whenever we write a set {ii, . . . ,ik}, it is implicitly understood that the elements 
ii,. . . ,ik are all distinct.) Clearly, each is an even element in the Grassmann 
algebra, and in particular it commutes with all the other elements of the Grassmann 
algebra. 

The definition (14. 2 p can also be rewritten as 



= Ul-\A\)+J2d^di)TA = [\{l~\A\) + ddyA 



(4.5) 

i,jeA 

where di = d/dipi and di = d/dipi are the traditional anticommuting differential 
operators satisfying diipj = 6ij, diipj = 0, diipj = Sij, diipj = and the (anti-)Leibniz 
rule, while d = ^ di and d = J2 ^i- 

Let us observe that 

{tavjb if|An5| = l 
f^A^rs = { (4.6) 

lO if\Ar]B\>2 

as an immediate consequence of (14.11) [when An B = {k}, only the second term in 
(14. 2 p with i = k survives]. Note, finally, the obvious relations 

hmi/i^) = {l-\A\)rA (4.7) 

A— >oo A 
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and 

= {X-m-\A\)TA. (4.8) 

We are interested in the subalgebra of the Grassmann algebra that is generated 
by the elements as A ranges over all nonempty subsets of V, for an arbitrary 
fixed value of Ajf] The key to understanding this subalgebra is the following amazing 
identity: 

Lemma 4.1 Let A,BCV with AnB ^ 0. Then 

A'^f'i' = r^"" ^f\^^Bi=i ^^^^^ 




More generally, 



/i^vr = r"' ■ (4.10) 

%f\AnB\>2 



if\AnB\>2 



/it;'] ^f\AnB\ = l 



where X" is the weighted average 

A" ^ - ^ ' . (4.n, 



(1^1 




-l)A + (|i?| 


\-l)X' i\A\ 




l)X + i\B\ 


\~1)X' 




1^1 


+ 


\B\ 




2 




\AUB\ 




1 



First proof. The formula f l4.10p can be proven by a direct (but lengthy) calculation 
within the Grassmann algebra that makes explicit a sort of fermionic-bosonic cancella- 
tion. Details can be found in the first preprint version of this article (arXiv : 0706 . 1509vl); 
see especially footnote 10 there. □ 

Second proof. We are grateful to an anonymous referee for suggesting the following 
simple and elegant proof using the differential operators d and d: 
Since d"^ = = 0, we have 

{ddTA) {ddTB) = dd [TAddTB) = Od {tbOBta) , (4.12) 



SO that 



/a ^ fP = X{1 - \A\)TAddTB + A'(l - \B\)TBddTA 

+ XX'il - \A\){1 - \B\)rATB + dd {TAddTB) . (4.13) 



^ One can also consider the smaller subalgebras generated by the elements f'^'' as A ranges over 
some collection S of subsets of V. 
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If \AnB\> 1, then taTb = and 

TAddTB = TBddTA = { ' ' (4.14) 

This proves (ICTl . □ 



j taub 


if 


\AnB\ 


= 1 




if 


\AnB\ 


> 2 



As a first consequence of Lemma [4.11 we have: 



(A) 



Corollary 4.2 Let A C V with \A\ > 2. Then the Grassmann element 
nilpotent of order 2, i.e. 

In particular, a product Hi^i fjt-^ vanishes whenever there are any repetitions 
among the Ai,...,Am. So we can henceforth concern ourselves with the case in 
which there are no repetitions; then E = {Ai,...,Am} is a set (as opposed to a 
multiset) and G = (V, E) is a hypergraph. 

By iterating Lemma 14.11 and using Proposition \2.2\ we easily obtain: 

Corollary 4.3 Let G = {V, E) be a connected hypergraph. Then 

T-r (\) I fv^ if G IS a hypertree 

UfA = { (4-15) 
A^E I if G is not a hypertree 



More generally, 



T-r (X.) \ fv*^ if G is a hypertree 

n /a = { (4.16) 
Ai^E I if G is not a hypertree 

where is the weighted average 

E(I^|-i)Aa E(I^I-i)Aa 

Age Age 

We are now ready to consider the subalgebra of the Grassmann algebra that is 
generated by the elements fA^ as A ranges over all nonempty subsets of V. Recall first 
that a partition of is a collection C = {G^} of disjoint nonempty subsets G^ O V 
that together cover V. We denote by Il{V) the set of partitions of V. If V has 
cardinality n, then n(V") has cardinality B{n), the n-th Bell number [37, pp. 33-34]. 
We remark that B{n) grows asymptotically roughly like n\ [38, Sections 6.1-6.3]. 

The following corollary specifies the most general product of factors f^\ Of 
course, there is no need to consider sets A of cardinality 1, since f^^^ = 1. 
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Corollary 4.4 Let E be a collection (possibly empty) of subsets ofV, each of cardi- 
nality > 2. 

(a) If the hypergraph G = {V,E) is a hyperforest, and {C^} is the partition of V 
induced by the decomposition of G into connected components, then Y\ f^^ = 

Y[fc^- More generally, Yl J'a^^ ~ Ylf'c'''^! where is the weighted average 
taken over the hyperedges contained in G^. 

(b) If the hypergraph G = {V,E) is not a hyperforest, then Yl fA^ = 0, and more 

Age 

generally Yl fj^"^^ — 0- 

A€E 

Proof. It suffices to apply Corollary 14 . 31 separately in each set G^, where {G^} is the 
partition of V induced by the decomposition of G into connected components. □ 

It follows from Corollary 14.41 that any polynomial (or power series) in the {f^^} 
can be written as a linear combination of the quantities f^'^^ = Yl^ flj^ for partitions 
C = {G^} e liiV). 

Using the foregoing results, we can simplify the Boltzmann weight associated to 
a Hamiltonian of the form 

n = - Y, ^Af^A^ ■ (4-18) 

Corollary 4.5 Let G = (y,E) be a hypergraph (that is, E is a collection of subsets 
ofV, each of cardinality >2). Then 

exp = (n-.^lriC). (4-19) 

\a&e J F e r{G) \AeF / Q=l 

where the sum runs over spanning hyperforests F in G with components Fi, . . . , Fi, 
and V{Fa) is the vertex set of the hypertree F^. More generally, 

expfj:^^') = E (U^^tlA-i'lr ("-20) 

\AeE / F&T{G) \AeF / a=l 

F = {Fu...,Fi) 

where is the weighted average i4-l'T\) taken over the hyperedges contained in the 
hypertree F^- 
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Proof. Since the f^^^ are nilpotent of order 2 and commuting, we have 

exp(5^u;^/?-M = ll{l + WAf['^^) (4.21a) 

XAeE ) A&E 



E'(1E \A&E' 



Using now Corollary 14.41 we see that the contribution is nonzero only when (V, E') is 
a hyperforest, and we obtain (HTT^ / fH^ . □ 



In a separate paper [39] we shall study in more detail the Grassmann subalgebra 
that is generated by the elements f^^^ as A ranges over all nonempty subsets of 
V . In the present section we have seen that any element of this subalgebra can 
be written as a linear combination of the quantities Z^^"* = f^^ for partitions 

C = {C^} G n(V^). It turns out that the quantities f^^^ are linearly dependent (i.e., 
an overcomplete set) as soon as \V\ > 4. We shall show [39], in fact, that a vector- 
space basis for the subalgebra in question is given by the quantities f^^^ as C ranges 
over all non-crossing partitions of V (relative to any fixed total ordering of V). It 
follows that the vector-space dimension of this subalgebra is given by the Catalan 
number Cn = where n = \V\. This is vastly smaller than the Bell number 

B{n), which is the dimension that the subspace would have if the f^'^^ were linearly 
independent. [Indeed, one can see immediately that the {fc'^^} must be linearly 
dependent for all sufficiently large n, simply because the entire Grassmann algebra 
has dimension only 4" <^ B{n).] It also turns out [39] that all the relations among 
the {/(^'^^} are generated (as an ideal) by the elementary relations Rabcd = 0, where 

-rtabcd — {a,b,c,d} ■> {b,c,d} ■> {a,c,d} ■> {a,b,d} ■> {a,b,c} 

, f{A) AX) , A\) A\) , .(A) .(A) . . r,r,^ 

J {aM-l {c,d} J{a,c}J{b,d} J{a,d}J{b,c} \^-^^) 

and a, 6, c, d are distinct vertices. 



5 Grassmann integrals for counting spanning hy- 
perforests 

For any subset A <ZV and any vector t = {ti)i^v of vertex weights, let us define 
the integration measure 

VaA^^.^) ■■= lld^.diJie'^^^'^'. (5.1) 
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Our principal goal in this section is to provide a combinatorial interpretation, in terms 
of spanning hyperforests, for the general Grassmann integral ( "partition function" ) 



Vi^ip, ip) exp 



Vv,t{ip,'(p) exp 



(5.2a) 



(5.2b) 



where G = {V, E) is an arbitrary hypergraph (that is, E is an arbitrary collection 
of subsets of V, each of cardinality > 2) and the {wa}a&e are arbitrary hyperedge 
weights. We also handle the slight generalization in which a separate parameter 
is used for each hyperedge A. 

Our basic results are valid for an arbitrary vector t = (ti)jgy of "mass terms". 
However, as we shall see, the formulae simplify notably if we specialize to the case in 
which ti = X for alH G ^. This is not an accident, as it corresponds to the case in 
which the action is 0SP(l|2)-invariant (see Section [7]). 

We begin with some formulae that allow us to integrate over the pairs of variables 
ipi, ipi one at a time: 



Lemma 5.1 Let A ^ V and i eV . Then: 

TA^i ifieA 
Uta ifi^A 



(b) I dtPid^iC'^^'"^^ ff^ 



A'li + iU - X)rA^, ifieA 



HJa 



ifi ^ A 



Proof, (a) is obvious, as is (b) when i ^ A. To prove (b) when z G A, we write 



Ia^ = ^(l-I^D^A + X^^Axj - 5Z '^j'^kTA^UM 



(5.3) 



J 6 A 



j, fc G A 



and integrate with respect to dipi dipi e''^^^'-^\ We obtain 

A(l - \A\)TA^i + tiTA^i + ^ TA-^{i,j} - i'ji'kTA 



(5.4) 



j, fc g A \ i 



(in the last term we must have j,k ^ i by parity), which equals + {U — X)ta-. 
as claimed. □ 



Applying Lemma [5. II repeatedly for i lying in an arbitrary set B (IV, we obtain: 
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Corollary 5.2 Let A,B CV. Then 



(A) 



n H fAiB + (E (^^-^) 

ieB^A / L iGBnA 



TA-^B 



In particular, for B = A we have 



A + Y.{t, - A) . 



(5.5) 



(5.6) 



i&A 



Proof. The factors ti for i G B^A follow trivially from the second line of Lemma ISTlT b). 
For the rest, we proceed by induction on the cardinality oi B r\ A. If |i? fl y4| = 0, 
the result is trivial. So assume that the result holds for a given set 5, and consider 
B' = B U {j} with j ^ A \ B. Using Lemma F5.1( a.b) we have 



\eBr\A ^ 

fflB' + ( E (^^ - ^)) 



(AxB)x{i} 



(5.7a) 
(5.7b) 



as claimed. □ 

Applying (15.61) once for each factor Ca, we have: 
Corollary 5.3 Let {Ca} be a partition ofV. Then 

vv,t{^, ^) n fc:^ = n + E - 



(5.8) 



The partition function (15.21) can now be computed immediately by combining 
Corollaries 14.51 and 15. 3[ We obtain the main result of this section: 

Theorem 5.4 Let G = (V, E) be a hypergraph, and let {waIage be hyperedge weights. 
Then 



T>{i{j, ip) exp 



Age 



F = {Fi,...,Fe) 



(5.9) 
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where the sum runs over spanning hyperforests F in G with components Fi, . . . , Fi, 
and V{Fa) is the vertex set of the hypertree Fa. In particular, if Xa takes the same 
value for all A, we have 



V{ip,'i/j) exp 



tiipiipi + ^ WAf^ 



AeE 



■ E n n h+ E c.-^) 

F G :F(G) \AeF / a=l \ ieV{F^) 
F=(Fu...,Fe) 

(5.10) 



Proof. We apply 05.81) . where (according to Corollary 14. 51) is the weighted average 
f l4.17p taken over the hyperedges contained in the hypertree Fa- Then 

+ (ti — Xa) 

ieV{F^) 



J2 t,-Xai\ViFa)\-l) 
ieViFc) 



ieV{Fc) AeEiFc) 



(5.11a) 
(5.11b) 



□ 



If we specialize (15.101) to tj = A for all vertices i, we obtain: 

Corollary 5.5 LetG = (V, E) he a hypergraph, and let {wA}AeE be hyperedge weights. 
Then 



V{i/j,ip) exp 



X^ilJiipi + ^WAf, 



(A) 
A 



Age 



FeJ^(G)\AeF J 



E n 

FGJf(G) VagF 



Wa 

Al^l- 



(5.12b) 



where the sum runs over spanning hyperforests F in G, and k{F) is the number of 
connected components of F. 

This is the generating function of unrooted spanning hyperforests, with a weight wa 
for each hyperedge A and a weight A for each connected component. Note that the 
second equality in f l5.12p uses Proposition 12. 1[ 

If, on the other hand, we specialize (15.101) to A = 0, we obtain: 

Corollary 5.6 Let G = (V, E) be a hypergraph, and let {wa}a€E be hyperedge weights. 
Then 



Viip^ip) exp 



AeE 



E n^n E«... 

F=iFi,...,Fe) 

(5.13) 
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where the sum runs over spanning hyperforests F in G with components Fi, . . . , Fi, 
and V{Fa) is the vertex set of the hypertree F^. 

This is the generating function of rooted spanning hyperforests, with a weight wa for 
each hyperedge A and a weight tj for each root i. 

Finally, returning to the case in which ti = X for all i, we can obtain a formula 
more general than fl5.12p in which the left-hand side contains an additional factor 
fc"^ = Yl fa^^ where C = {Cy} is an arbitrary family of disjoint nonempty subsets of 

V. Indeed, it suffices to differentiate fl5.12p with respect to all the weights ifc^,!^ We 
obtain: 

Corollary 5.7 Let G = (y,E) be a hypergraph, let {wa}a<=e be hyperedge weights, 
and let C = {C-y} be a family of disjoint nonempty subsets ofV. Then 

i A&E 



= E (n 

F&T{G;C)\AeF 

where J-'{G; C) denotes the set of spanning hyperforests in G that do not contain 
any of the {C^} as hyperedges and that remain hyperforests (i.e., acyclic) when the 
hyperedges {G^} are adjoined. 

Indeed, to deduce Corollary 15 . 71 from Corollary I5.5l by differentiation, it suffices to ob- 
serve that, by Proposition [2111 the number of connected components in the hyperforest 
obtained from F by adjoining the hyperedges {G^} is precisely k{F) — J2i\^-y\ ~ 

7 

For instance, if H /c^^ consists of a single factor fc, then {G}) consists of 

7 

the spanning hyperforests in which all the vertices of the set G belong to different 

components. Similarly, if H /c"^ consists of two factors fcifc2 with Ci fl C2 = 0, 

7 

then J-'{G; {Ci, C2}) consists of the spanning hyperforests in which each component 
contains at most one vertex from Gi and at most one vertex from C2. The conditions 
get somewhat more complicated when there are three or more sets G^. 

It is possible to obtain an analogous extension of Theorem 15. 41 by the same method, 
but the weights get somewhat complicated, precisely because we lose the opportunity 
of using Proposition 12.11 in a simple way. 

^ If the sets C-y do not happen to belong to the hyperedge set E, it suffices to adjoin them to 
E and give them weight wc^ = 0. Indeed, there is no loss of generality in assuming that G is 
the complete hypergraph on the vertex set i.e. that every subset of V of cardinality > 2 is a 
hyperedge. 



|V{^P,^P) m/SM exp 



Wa 



(5.14) 
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Equations fl5.9l) - fl5.13l) are the hypergraph generahzation of fll.5l) - fll.8p . respec- 
tively. To see this, let G = (V, E) be an ordinary graph, so that each edge e E E 
is simply an unordered pair {i,j} of distinct vertices i,j G V, to which there is 
associated an edge weight Wij = wji. Then by definition (14.21) we have 



so that if we take Ajj 



-Uij we have 



{i,j}&E i,j€V {i,j}eE 



(5.15) 
(5.16) 



(5.17) 



where the (weighted) Laplacian matrix for the graph G is defined as 

'Wij iii^ j 
J2 Wik if i = j 

Then f l5.9p - fl5.13p become precisely f ll.5p -f |T78l) . respectively. 

More generally, consider the case in which G = {V, E) is a /c-uniform hypergraph 
(an ordinary graph corresponds to the case k = 2). Let Wi^^,,,^i^ (assumed com- 
pletely symmetric in the indices ii, . . . ,ik) be the weight associated to the hyperedge 
{ii, . . . , ik} when ii, . . . ,ik are all distinct, and let Wi^^,,,^i^ = when at least two in- 
dices are equal. Define the (weighted) Laplacian tensor (a rank-fc symmetric tensor) 
by 



—w 



ti,...,ik 



Then we have 



if 4 are all different 

if ir = is {t ^ s) and the others are all different 

otherwise 

(5.18) 



E 



ii,...,ik£V 



{k-2)\ 



(5.19) 

so that the "action" is given by f l5.19p plus the "mass term" XY^-ipiipi. Combining 
Corollary 15. 51 with fl5.19p . we obtain a formula for the generating function of spanning 
hyperforests in a fc-uniform hypergraph: 



XHF) 



1 v,*fc 



il,...,jftGV 



ik-2)\ 



A -r 



(5.20) 
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Let us remark that while the Laplacian matrix fl5.17l) for an ordinary graph has 
vanishing row and column sums (i.e., Ylj-^ij — 0); Laplacian tensor flS.lSp for 
a hypergraph satisfies Li-^^,,,^i^ = when ii, . . . ,ik-i are all distinct, but not in 
general otherwise. 

For an application to counting spanning hyperforests in the complete /c-uniform 
hypergraph, see [40]. 

Remark. One can also generalize (15. 2p to allow products /^^^ = Yl fc"^ the 

a 

exponential (i.e., in the action) in place of the single operators with correspond- 
ing coefficients wc- These generalized integrals likewise lead to polynomials in the 
variables {wc} such that the union of the families Cj arising in any given monomial 
is the set of hyperedges of a hyperforest. However, the simultaneous presence of 
certain sets of hyperedges in the hyperforest now gets extra weights. We hope to 
discuss these extensions elsewhere. This generalized model is conceptually important 
because, when ti = X for all i, it corresponds to the most general 0SP(l|2)-invariant 
action (see Section [Tj). 



6 Extension to correlation functions 

In the preceding section we saw how the partition function (15. 2p of a particular 
class of fermionic theories can be given a combinatorial interpretation as an expan- 
sion over spanning hyperforests in a hypergraph. In this section we will extend this 
result to give a combinatorial interpretation for a class of Grassmann integrals that 
correspond to (unnormalized) correlation functions in this same fermionic theory; we 
will obtain a sum over partially rooted spanning hyperforests satisfying particular 
connection conditions. 

Given ordered fc-tuples of vertices / = (ii, Z2, • • • , ik) & and J = (ji, j2, • • • , jk) £ 
V'', let us define the operator 

Oi,j := ^hipn-'-At-^jk . (6-1) 

which is an even element of the Grassmann algebra. Of course, the ii, i2, ■ ■ ■ ,ik must 
be all distinct, as must the ji, j2, • • • ,jk, or else we will have Oi^j = 0. We shall 
therefore assume henceforth that I, J ^ V^, where is the set of ordered fc-tuples 
of distinct vertices in V. Note, however, that there can be overlaps between the 
sets {ii,i2, • • • , h} and {ji, j2, • • • ,jk}- Note finally that Ojj is antisymmetric under 
permutations of the sequences / and J, in the sense that 

OioajoT = sgn((T) sgn(r) (6.2) 

for any permutations cr, r of {1, . . . , k}. 
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Our goal in this section is to provide a combinatorial interpretation, in terms of 
partially rooted spanning hyperforests satisfying suitable connection conditions, for 
the general Grassmann integral ( "unnormalized correlation function" ) 



[Oj^j] = Z{Ojj) 



Vi^j,^) Ojj exp 



^ tilpilpi + ^ WAf, 



AeE 



Vv,t{ip,ip) Oij exp 



(A) 
A 



.A€E 



(A) 



i€A 

1 ifk = l 







ifk>2 



(6.3a) 



(6.3b) 



The principal tool is the following generalization of (15.61) : 

Lemma 6.1 Let A C V, and let I = {ii,i2, ■ ■ ■ ,ik) ^ (ind J = (ji, j2, • • • ,jk) ^ 
A^. Then 



(6.4) 



Proof. The case k = is just (15.61) . To handle k = 1, recall that 

m,'^A\{l,m} 



J A 



\{1-\A\)ta + J^TA-., 



(6.5) 



Now multiply f^^^ by with j G A, and integrate with respect to X'A,t(V^, '^)- If 
i = j, then the only nonzero contribution comes from the term i = i in the single sum, 
and ipiipiTA-^i = ta, so the integral is 1. If z 7^ j, then the only nonzero contribution 
comes from the term i = j, m = i in the double sum, and (^/'i^/'j)(— ^/'j^/'j)rA\{ij} = ta, 
so the integral is again 1. 

Finally, if |/| = |J| = A; > 2, then every monomial in Oj^jf^^ has degree > 
2\A\ - 2 + 2k>2\A\, so Oi^jf^J"^ = 0. □ 

Of course, it goes without saying that if m{ip, ip) is a monomial of degree k in the 
variables ipi {i G A) and degree k' in the variables ipi (z G A), and k is not equal to 
k', then JVA,t{^,tlj)m{ij,^/j)f['^ = 0. 



Now go back to the general case / = (ii, ^2, • • • , ik) £ and J = (ji, j2, ■ ■ ■ ,jk) G 
V^, let C = {Ca}a=i a partition of V, and consider the integral 



T(/,J;C) 



D^,t(^,^)07,jn/cj 

a=l 



(6.6) 
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The integral factorizes over the sets Ca of the partition, and it vanishes unless |/ fl 
Cq\ = I JHCqI for all a; here lP\Ca denotes the subsequence of / consisting of those 
elements that lie in Ca, kept in their original order, and | / fl Cq,| denotes the length 
of that subsequence (and likewise for JflCa). So let us decompose the operator Oi^j 
as 

inCcJnc^ , (6.7) 

a=l 

where a{I,J;C) G {±1} is a sign coming from the reordering of the fields in the 
product. Applying Lemma [6. II once for each factor Ca, we see that the integral (16.61) 
is nonvanishing only if |/ fl Cq,| = | J fl Ca| < 1 for all a: that is, each set Ca must 
contain either one element from I and one element from J (possibly the same element) 
or else no element from I or J. Let us call the partition C properly matched for (J, J) 
when this is the case. (Note that this requires in particular that m > k.) Note 
also that for properly matched partitions C we can express the combinatorial sign 
a{I, J; C) in a simpler way: it is the sign of the unique permutation vr of {1, . . . ,k} 
such that ir and j,r(r) lie in the same set Ca for each r {1 < r < k). (Note in particular 
that when {ii, 12, ■ ■ ■ , ik} H {ji,j2, ■ ■ ■ ,jk} = S 0, the pairing vr has to match the 
repeated elements [i.e., v = jn(r) whenever v € S], since a vertex cannot belong 
simultaneously to two distinct blocks Ca and Cp.) We then deduce immediately from 
Lemma 16.11 the following generalization of Corollary 15.31 

Corollary 6.2 Let I, J & and let C = {Ca} be a partition of V. Then 

a 

sgn(7r) Y[ ( A + X] (^i ^ A) ) if C is properly matched for (/, J) 

a:\InCc,\=0^ i&Cc, ^ (6.8) 

otherwise 

where vr is the permutation 0/ {1, . . . , fc} such that ir and ^^^(r) He in the same set Ca 
for each r. 

We can now compute the integral (16. 3p by combining Corollaries 14.51 and 16. 2[ If 
C = (V, E) is a hypergraph and C is a spanning subhypergraph of G, let us say that 
C is properly matched for (/, J) [we denote this by C ~ (J, J)] in case the partition of 
V induced by the decomposition of C into connected components is properly matched 
for (/, J). We then obtain the main result of this section: 
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Theorem 6.3 Let G = (y,E) be a hypergraph, let {wA}AeE be hyperedge weights, 
and let I,JeV^. Then 



V{ij,^p) Oij exp 



p(A) 
I A 



i AeE 

= ^Snircj^j.,^) ( n ^a) n + E - ^)) ' (6.9) 

F ~ (7, J) 
F= (Fi,...,F,) 

where the sum runs over spanning hyperforests F in G, with components Fi, . . . , Fi, 
that are properly matched for (/, J), and V{Fa) is the vertex set of the hypertree F^; 
here t^i,j-f is the permutation of {l^...,k} such that v o-nd ^^^(r) lie in the same 
component Fa for each r. 

If we specialize fl6.9p to tj = A for all vertices i, we obtain: 

Corollary 6.4 Let G = (y,E) be a hypergraph, let {wA}AeE be hyperedge weights, 
and let I,JeV^. Then 



V{i,,i,) Oi,j exp 



(A) 
A 



x^i^ii^i + ^WAf, 

i AeE 

e r(G) \AeF J 



F e TiG) 
F~(/,J) 



Vl-fc 



E 

F e ^(G) 

F~(/,J) 



n 

KAeF 



Wa 



X 



A\-l 



(6.10a) 



(6.10b) 



where the sum runs over spanning hyperforests F in G that are properly matched 
for (/, J), and k{F) is the number of connected components of F; here t^i^j-f is the 
permutation of {1, ... ,k} such that i^ and j,r(r) in the same component of F for 
each r. 

This is the generating function of spanning hyperforests that are rooted at the vertices 
in J, J and are otherwise unrooted, with a weight wa for each hyperedge A and a 
weight A for each unrooted connected component. 

If, on the other hand, we specialize (16.91) to A = 0, we obtain: 
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Corollary 6.5 Let G = {V,E) be a hypergraph, let {wa}a€E be hyperedge weights, 



and let I,J eV^. Th en 



V{iP,i,) Oi,j exp 



i AeE 

sgn(7r,,j;p) ( n n ( E ' ^^-^^^ 



Fejc-(G) V^eF / a: \inV{F„)\=o \eV{F„) 

F ~ (/, J) 



where the sum runs over spanning hyperforests F in G, with components Fi, . . . , Fi, 
that are properly matched for (/, J), and V{Fa) is the vertex set of the hypertree Fa; 
here t^i,j;f is the permutation of such that v (^nd j^^i^j.) lie in the same 

component F^ for each r. 

This is the generating function of rooted spanning hyperforests, with a weight wa for 
each hyperedge A and a weight ti for each root i other than those in the sets /, J. 

Let us conclude by making some remarks about the normalized correlation func- 
tion {Ojj) obtained by dividing fl6.3l) by fl5.2l) . For simphcity, let us consider only 
the two-point function (ipitpj). We have 



7., (a+ (tk-X))) , (6.12) 



\ ^ fcer(i) 

where the expectation value on the right-hand side is taken with respect to the 
"probability distribution"^ on spanning hyperforests of G in which the hyperforest 
F = (Fi, . . . , Ff) gets weight 

\AeF J a=l keV{Fc) 

•jij denotes the indicator function 

{1 if i and j belong to the same component of F 
it not ^'■'"^ 

and r{i) denotes the vertex set of the component of F containing i. The factor 

( A + J2 i^k — X)] in (16.121) arises from the fact that in fIS.lOp each component 
^ fcer(i) ^ 



We write "probability distribution" in quotation marks because the "probabilities" will in 
general be complex. They will be true probabilities (i.e., real numbers between and 1) if the 
hyperedge weights wa are nonnegative real numbers. 
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gets a weight A + '^kGr{i)i^k ~ '^)' while in (16. 9p only those components other than 
the one containing i and j get such a weight. So in general the correlation function 
(iJii/jj) is not simply equal to (or proportional to) the connection probability (7ij). 
However, in the special case of Corollaries 15.51 and 16.41 — namely, all ti = A, so that 
we get unrooted spanning hyperforests with a "flat" weight A for each component — 
then we have the simple identity 

(V^,^,) = A-i(7^,) . (6.15) 

Combinatorial identities generalizing (16.151) . and their relation to the Ward identities 
arising from the 0SP(1|2) supersymmetry, will be discussed elsewhere [31]. 



7 The role of 0SP(1|2) symmetry 

In [12] we have shown how the fermionic theory (ll.7p emerges naturally from the 
expansion of a theory with bosons and fermions taking values in the unit supersphere 
in M^l^, when the action is quadratic and invariant under rotations in 0SP(1|2). Here 
we would like to discuss this fact in greater detail, and extend it to the hypergraph 
fermionic model (I5.12p . 

We begin by introducing, at each vertex i E V, a superfield iij := {ai,ipi,ijji) 
consisting of a bosonic (i.e., real) variable cxj and a pair of Grassmann variables 
tjji,i'i- We equip the "superspace" M^'^ with the scalar product 

n, ■ Uj := aiCj + Xi^pi^pj - V'i^j) , (7.1) 

where A 7^ is an arbitrary real parameter. 

The infinitesimal rotations in M^'^ that leave invariant the scalar product (17.10 
form the Lie superalgebra osp(l|2) [41-43]. This algebra is generated by two types of 
transformations: Firstly, we have the elements of the sp(2) subalgebra, which act on 
the field as = rij + Sui with 

Sai = (7.2a) 
Sipi = -atpi + ^fiJi (7.2b) 
Sipi = +ai)i + (3i)i (7.2c) 

where ce,f3,^ are bosonic (Grassmann-even) global parameters; it is easily checked 
that these transformations leave (17. ip invariant. Secondly, we have the transforma- 
tions parametrized by fermionic (Grassmann-odd) global parameters e, e: 

Sa, = -xy^(e^^ + ^p^e) (7.3a) 
6iJi = X-^^\ai (7.3b) 
6iJi = A^^/^ecTi (7.3c) 
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(Here an overall factor A~^/^ has been extracted from the fermionic parameters for 
future convenience.) To check that these transformations leave (17.11) invariant, we 
compute 

6{ni ■ Uj) = {6(Ji)aj + (Ti{6aj) + A [(^V'OV'i + i'ii^'^j) - - '^ii^i'j)] 

(7.4a) 

+ A^/^ [eipjai + ipieaj - eipjai - tpiiaj] (7.4b) 
= . (7.4c) 

In terms of the differential operators di = d/dipi and di = d/dipi, the transforma- 
tions (17.21) can be represented by the generators 



^0 = X^(V^i<9i - ipidi) (7.5a) 
X+ = ^^,9, (7.5b) 

X_ = Y.'^A (T.5c) 

iev 

corresponding to the parameters a, f5, 7, respectively, while the transformations (17.31) 
can be represented by the generators 

corresponding to the parameters e, e, respectively. (With respect to the notations 
of [43] we have X± = L=p, Xq = — 2Lo and Q± = ^2iR^.) These transformations 
satisfy the commutation/anticommutation relations 

[Xo,X±] = ±2X± [X+,X_]=Xo (7.7a) 

{Q±,Q±} = ±2X± {g+,g_} = Xo (7.7b) 

[Xo,Q±] = ±g± [x±,g±] = [x±,Q^] = -g± (7.7c) 

Note in particular that X± = Q\ and Xg = Q+Q^ + g_g+. It follows that any 
element of the Grassmann algebra that is annihilated by Q± is also annihilated by 
the entire osp(l|2) algebra. 
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Now let us consider a cr-model in which the superfields rij are constrained to he 
on the unit supersphere in i.e. to satisfy the constraint 

Ui-m = a^ + 2Xtljiiji = 1. (7.8) 

We can solve this constraint by writing 

ai = ±(1 - 2A^,^,)^/2 _ ±(i_A^,^,), (7.9) 

exploiting the fact that ipf = ipf = 0. Let us henceforth take only the + sign in (17. 9p . 
neglecting the other solution (the role played by these neglected Ising variables will 
be considered in more detail elsewhere [30]), so that 

ai = 1- X^i^, . (7.10) 

We then have a purely fermionic model with variables ipyip in which the sp(2) trans- 
formations continue to act as in 07.21) while the fermionic transformations act via the 
"hidden" supersymmetry 

S^|Ji = A-^/2^(l-A^,^,) (7.11a) 
5^, = A-i/2^(l-A^,^,) (7.11b) 

All of these transformations leave invariant the scalar product 

ni -Uj = 1 - X{ipi - - ipj) + X'^ipiipiipjipj . (7.12) 

The generators Q± are now defined as 

Q+ = X-'/'J2il-X^P,^,)^, = X-'/^d - X'/'J^'^.m (7.13a) 

iev i&v 

Q_ = X-y^J2i'^-X^p,^p,)^, = X-'/^d - X'/'J2'^,^PA (7.13b) 

iGV i&V 

where we recall the notations d = 'Y^ di and d = ^i- 

i£V i&V 

Let us now show that the polynomials f^^^ defined as in (14.51) are 0SP(1|2)- 
invariant, i.e. are annihilated by all elements of the osp(l|2) algebra. As noted previ- 
ously, it suffices to show that the f^^^ are annihilated by Q±. Applying the definitions 
f l7.13p . we have 

Q-TA = X-^'^dTA (7.14) 

and hence 

Q+Q-TA = X-'ddTA - \A\ta , (7.15) 

so that 

/f = A(l + g+Q_)r^. (7.16) 
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The next step is to compute Q+f^^: since 

g+(l + Q+Q_) = Q+ + QIQ- = Q+ + X+Q. 

= Q+ + [x+, g_] + g_x+ = Q+-Q+ + g_x+ = g_x+ (7.17) 



by the relations ( ]7.7bl) / ( !7.7cp . while it is obvious that X+r^ = 0, we conclude that 



Q+f^^ = 0, i.e. is invariant under the transformation g+. A similar calculation 
of course works for g_l"1 

In fact, the 0SP(l|2)-invariance of can be proven in a simpler way by writing 
explicitly in terms of the scalar products iij ■ rij for i,j G A. Note first that 

43} " -X'lpi'ipi'ipj'ipj + {-ipi - - -ipj) (7.18a) 

= l(l-n,-n,) (7.18b) 



A 
2A 



(7.18c) 



By Corollary 14.31 we obtain 



= (1 - ■ ^i2) (1 - ■ ^h) ■ • ■ (1 - ni^_, ■ Hi J (7.19 



1 



a 



(2A) 



— (n,, - riij (n,2 - n^g) ■ ■ ■ (n,^_^ - n^J . (7.19b) 



Note the striking fact that the right-hand side of (17.191) is invariant under all per- 
mutations oi ii, . . . ,ik, though this fact is not obvious from the formulae given, and 
is indeed false for vectors in Euclidean space with N ^ —1. Moreover, the 
path ii, . . . ,ik that is implicit in the right-hand side of fl7.19p could be replaced by 
any tree on the vertex set {ii, . . . ,ik}, and the result would again be the same (by 
Corollary US]). 

It follows from (17. 18p / (I7.19P that the subalgebra generated by the scalar products 
rij ■ rij for i,j & V is identical with the subalgebra generated by the for A (1 V, for 
any A 7^ 0. Therefore, the most general 0SP(l|2)-symmetric Hamiltonian depending 
on the {rijjjgy is precisely the one discussed in the Remark at the end of Section [5l 
namely in which the action contains all possible products f^^^ = Ylfc'^y where {Ca} 

a 

is a partition of V. 
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We are grateful to an anonymous referee for suggesting this proof. An alternate proof that 
Q±fA^ — 0; based on direct calculation using the definition glU) of f'^^\ can be found in the first 



preprint version of this article (arXiv : 0706 . 1509vl): see equations (7.8)-(7.11) there 
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In Appendix El we will prove a beautiful alternative formula for f^^^^^ 



fu\ ; = TTTrrdetM (7.20) 

where M is the k x k matrix of scalar products Mrs = ■ rij^. In this formula, 
unlike f l7.19p . the symmetry under all permutations of ii,...,ik is manifest. We 
remark that the determinant of a matrix of inner products is commonly called a 
Gram determinant [44, p. 110]. 

Finally, we need to consider the behavior of the integration measure in fl5.2p . 
namely 

1^v,t{4^,^P) = (7-21) 
iev 

under the supersymmetry (17.111) . In general this measure is not invariant under 
( 17. lip , but in the special case ti = X for all i, it is invariant, in the sense that 

Vv,xiiJ,^p)SF{i;,iJ) = (7.22) 

for any function F{il),%l>). Indeed, T>v,x{ip,ip) is invariant more generally under local 
supersymmetry transformations in which separate generators ej,ej are used at each 
vertex i. To see this, let us focus on one site i and write -F(V', ip) = a + bipi + cipi + 
dipiipi where a,b,c,d are polynomials in the (which may contain both 

Grassmann-even and Grassmann-odd terms). Then 

6F = X^^'^ [beiCTi + ciiai + d [eiaitpi + tpieiai)] (7.23a) 



1 



0"i 



^1/2 



hei + cei + d(eiVi + ^i^) 



(7.23b) 



Since ai = e~^'^^'^\ this cancels the factor e^^^'^' from the measure (since ti = X) 
and the integral over dipidipi is zero (because there are no ipiipi monomials). Thus, 
the measure Vy^ti'^,'^) is invariant under the local supersymmetry at site i when- 
ever ti = X. If this occurs for all i, then the measure is invariant under the global 
supersymmetry ( 17. 3p . 

The 0SP(l|2)-invariance of T>v,x{ip-,ip) can be seen more easily by writing the 
manifestly invariant combination 

(5(n- - 1) drii = 6{a'^ + 2Xi)i^i - 1) dai dipi dipi (7.24a) 
= e^^^^^ 6{ai - (1 - XiJiipi)) dai dipi d^i , (7.24b) 

where the factor e^"^^"^^ comes from the inverse Jacobian. Integrating out cxj from 
(l7.24bD . we obtain e^'^''^' d^i d^i. 
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As a consequence of (]7.19bl) and (]7.24b|) . the generating function (15.201) for span- 
ning hyperforests in a fc-uniform hypergraph can be rewritten as 



— y 



— 1) cin,- X 



exp 



^ ^ ii,...,%k&V ^ 



. (7.25) 



In the special case k = 2, this result appears in [12]. 



8 Conclusions 

In this paper we have applied techniques of Grassmann algebra, first used in [12] 
to obtain the generating function of spanning forests in a graph — generalizing Kirch- 
hoff 's matrix-tree theorem — to a wider class of models associated to hypergraphs. 
The key role in our analysis is played by a set of simple algebraic rules (Lemma 14.11) 
that express, in a certain sense, the fermionic-bosonic cancellation associated to the 
underlying 0SP(1|2) supersymmetry. This algebraic approach allows for notably 
simplified proofs and for strong generalizations. 

In particular, we are able to obtain combinatorial interpretations in terms of 
spanning hyperforests for the partition function (Section |5]) and correlation functions 
(Section [6]) of a fairly general class of fermionic models living on hypergraphs. In 
that subset of the results where the 0SP(1|2) supersymmetry is preserved, the com- 
binatorial weights of the hyperforest configurations become (perhaps not surprisingly) 
notably simpler. Among other things, we obtain the generating function of unrooted 
spanning hyperforests on a weighted hypergraph (together with a family of relevant 
combinatorial observables) as an 0SP(l|2)-invariant fermionic integral. 

Finally, in Appendix [B] we present a graphical formalism for proving both the 
classical matrix-tree theorem and numerous extensions thereof, which can serve as an 
alternative to the algebraic approach used in the main body of this paper and which 
we hope will have further applications. 

In a follow-up paper [39] we shall study in more detail the Grassmann subalgebra 
that is generated by the elements as A ranges over all nonempty subsets of V. 

It is also natural to ask about extensions of this work in which combinatorial inter- 
pretations are obtained for statistical-mechanical models with other supersymmetry 
groups. We are currently studying models with 0SP(l|2n) and OSP(2|2) supersym- 
metries and hope to report the results in the near future. 
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A A determinantal formula for f)^ 

The main purpose of this appendix is to prove the determinantal formula (17.201) 
for Along the way we will obtain a rather more general graphical representation 
of certain determinants. 

Let A = {aij)^j^i be a matrix whose elements belong to a commutative ring R. 
The determinant is defined as usual by 

n 

detA = ^ sgn(7r) JJai^(i) , (A.l) 

7rgn,i i=l 

where the sum runs over permutations vr of [n] := {1, . . . , n}, and 
sgn(7r) = (_i)#(even cycles of tt) ^^iq sign of the pcrmutatiou vr. 

We begin with a formula for the determinant of the sum of two matrices in terms 
of minors, which ought to be well known but apparently is not0: 

Lemma A.l Let A and B benxn matrices whose elements belong to a commutative 
ring R. Thei^ 

det{A + B) = e{I,J){detAij){detBicjc), (A.2) 

I,J C [n] 
\I\ = \J\ 

where e{I,J) = (— l)Sie/«+EjejJ is the sign of the permutation that takes IP into 
JJ^ (where the sets I,P, J, P are all written in increasing order). 

Proof. Using the definition of determinant and expanding the products, we have 
det{A + B) = sgn(7r) ^ JJ ai^^i) JJ . (A.3) 

Define now J = 7i[I]. Then we can interchange the order of summation: 

det{A + B) = Y Yl sgn(7r) JJ ai^(i) JJ . (A.4) 

I,J C[n] TT G n„ i^I f^&I" 

\I\ = \ J\ tt[I] = J 



This formula can be found in [45, pp. 162-163, Exercise 6] and [46, pp. 221-223]. It can also be 
found — albeit in an ugly notation that obscures what is going on — in [45, pp. 145-146 and 163- 
164] [47, pp. 31-33] [48, pp. 281-282]; and in an even more obscure notation in [49, p. 102, item 5]. 
We remark that an analogous formula holds (with the same proof) in which all three occurrences of 
determinant are replaced by permanent and the factor e(/, J) is omitted. 

The determinant of an empty matrix is of course defined to be 1. This makes sense in the 
present context even if the ring R lacks an identity element: the term I = J = contributes det_B 
to the sum (|A.2|) . while the term I = J = [n] contributes detA. 
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Suppose now that |/| = \ J\ = and let us write / = {ii, . . . , ik] and J = {ji, . . . 
where the elements are written in increasing order, and likewise = {ii, . . . ,in-k} 
and J = {mi, . . . ,m„_fc}. Let tt' G 11^ and tt" G Hn-k be the permutations defined 
so that 

T^{ia) = jp < ^ 7r'(a) = (3 (A.5a) 

ni^Ca) = mp < — > 7i"{a) = f3 (A. 5b) 

It is easy to see that sgn(7r) = sgn(7r') sgn(7r")e(/, J). The formula then follows by 
using twice again the definition of determinant. □ 



Corollary A. 2 Let A and B be n x n matrices whose elements belong to a commu- 
tative ring R. Then det{A + XB) is a polynomial in X of degree at most rank(_B), 
where "rank" here means determinantal rank (i.e. the order of the largest nonvanish- 
ing minor). 

Proof. This is an immediate consequence of the formula (1A.2I) . since all minors of 
B of size larger than its rank vanish by definition. □ 

Next recall the traditional graphical representation of the determinant: 

Lemma A. 3 Let C = {cij)'^j^i be a matrix whose elements belong to a commutative 
ring R. Then 

det(-C) = ^(_l)#(cyclesofG) -Q ^ (A g) 

G ijeE{G) 

where the sum runs over all permutation digraphs G on the vertex set {1,2, . . . ,n}, 
i.e. all directed graphs in which each connected component is a directed cycle (possibly 
of length 1 ). 

Proof. This is an immediate consequence of flA.ll) and the fact that (-1)#('=™'^ '=y'=i'=« °^ 

^ -j^ ^ # (cycles of 7r) ^ -|^^ # (odd cycles of tt) | | 

Now let a = ^^^d b = be a pair of vectors with elements in the ring 

R. The main result of this appendix is the following generalization of Lemma IA.3I 

Lemma A. 4 Let G = {cij)^^^^ be a matrix whose elements belong to a commutative 
ring R, and let a = (aj)"^]^ and h = be vectors with elements in R. Then 

det(abT-C) = det(-C) + H ^^^^ ' (A'^) 

G ij<^E{G) 

where the sum runs over all directed graphs G on the vertex set {1, 2, . . . , n} in which 
one connected component is a directed path (possibly of length 0, i.e. an isolated vertex) 
from source s{G) to sink t{G) and all the other connected components are directed 
cycles (possibly of length 1). 
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Proof. Introduce an indeterminate A and let us compute det(Aab"'" — C), working 
in the polynomial ring -R[A], by substituting Cij — Xaibj in place of Cij in (1A.6I) . The 
term of order A° is det(— C), which is given by (]A.6I) . In the term of order A^, one 
edge ij in G carries a factor — aj6j and the rest carry matrix elements of C. Setting 
G' = G \ ij , we see that G' has one less cycle than G [thereby cancelling the minus 
sign] and has a path running from source s{G') = j to sink t{G') = i. Dropping the 
prime gives (lA.TI) . Terms of order A^ and higher vanish by Corollary I A. 21 because ab""" 
has rank 1. □ 

Corollary A. 5 Let G = {cij)^j^i be a matrix whose elements belong to a commuta- 
tive ring-with-identity element R, and let E be the n x n matrix with all elements 1. 
Then 

det(E- C) = det(-C) + Yl c^j , (A.8) 

where the sum runs over all directed graphs G on the vertex set {1, 2, . . . , n} in which 
one connected component is a directed path (possibly of length 0, i.e. an isolated vertex) 
and all the other connected components are directed cycles (possibly of length 1). 

The following result is an immediate consequence of Lemma fA.SI and Corollary lA. 51 

Corollary A. 6 Let G = {cij)^j^^ be a matrix whose elements belong to a com- 
mutative ring-with-identity- element R and satisfy Ci^i2Ci2i-j, ■ ■ ■ Ci^_iik(^ikii = for all 
ii, . . . ,ik {k > 1), and let E be the n x n matrix with all elements 1. Then 

det{E -G) = J2 n ' (^-9) 

P ijGE{P) 

where the sum runs over all directed paths P on the vertex set {1, 2, . . . , n}. (There 
are n\ such contributions.) 

Proof. The hypotheses on G lead to the vanishing of all terms containing at least 
one cycle (including cycles of length 1). Therefore, the only remaining possibility is 
a single directed path. □ 

Let us now specialize Corollary IA.6I to the case in which the commutative ring R 
is the even subalgebra of our Grassmann algebra, and the matrix G is given by 

Cii = (A. 10) 

Cij = Cji = A/|2j for i ^ j (A. 11) 



36 



The hypothesis Cj^jjCjjig • • • Cjj,_^jj.Cj^j^ = is an immediate consequence of Corol- 
lary |l]3l Moreover, by equation ( ]7.18bl) we have {E — C)ij = n, ■ iij. In the expansion 
flA.QP we obtain n! terms, each of which is of the form A"~^ times nijeE(P) 
some directed path P on the vertex set {1,2, ... ,n}. But by Corollary 14.31 each such 
product equals /|^^ ^y, so this proves the determinantal formula (I7.20p for /^^^ 

B Graphical proof of some generalized matrix-tree 
theorems 

In this appendix we shall give a "graphical" proof of the classical matrix-tree 
theorem as well as a number of extensions thereof, by interpreting in a graphical way 
the terms of a formal Taylor expansion of an action belonging to the even subalgebra 
of a Grassmann algebra. (We require the action to belong to the even subalgebra in 
order to avoid ordering ambiguities when exponentiating a sum of terms.) Some of 
these extensions of the matrix-tree theorem are already set forth in the main body of 
this paper, where they are proven by an "algebraic" method based on Lemma SiT] and 
its corollaries. Other more exotic extensions are described here with an eye to future 
work; they could also be proven by suitable variants of the algebraic technique. 

Curiously enough, it turns out that the more general is the fact we want to prove, 
the easier is the proof; indeed, the most general facts ultimately become almost 
tautologies on the rules of Grassmann algebra and integration. The only extra feature 
of the most general facts is that the "zoo" of graphical combinatorial objects has to 
become wider (and wilder). 

So, in this exposition we shall start by describing the most general situation, and 
then show how, when special cases are chosen for the parameters in the action, a 
corresponding simplification occurs also in the combinatorial interpretation. 

B.l General result 

Consider a hypergraph G = {V, E) as defined in Section [21 i.e. V" is a finite set 
and E is a. set of subsets of V , each of cardinality at least 2, called hyperedges. As 
usual we introduce a pair of Grassmann generators for each i & V. We shall 

consider actions of the form 

5(^,V) = J]5a(^,^), (B.l) 

where 

5a (^,^) = w\Ta + "^WA^iTA-^i + ^ WA;ijA'4'j'TA~^{i,j} (B.2) 

i,j€A 
i 7^ J 
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and ta_ = YlieA'^i'^i- Please note that the form (1B.2P resembles the definition fl4.2p of 
the same monomials appear, but now each one is multiplied by an independent 
indeterminate. Thus, for each hyperedge A of cardinality k we have + 1 parameters: 
w\, {wA-i}ieA and {wA;i,j}(i^j)£A- [We have chosen, for future convenience, to write 
the last term in (14.21) as +ipiipj rather than —ipiiljj] 

Please note that, for \A\ > 2, all pairs of terms in Sa{iPi'4') have a vanishing 
product, because they contain at least 2(2|yl| — 2) = 4|yl| —4 fermions in a subalgebra 
(over A) that has only 21^41 distinct fermions. As a consequence, we have in this case 

exp[5^(V^,^)] = l + 5^(7/;,7/'). (B.3) 

On the other hand, if \A\ = 2 (say, A = {i,j}), we have two nonvanishing cross- 
terms: 

(wA-i'4'jipj) {wA;ji^ii^i) = WA-iWA-j^ii^ii'ji^j (B.4a) 

{wA;ijAi'j) i^A-Ji^ji'i) = -WA;ijWA;jiAA4'j^j (B.4b) 

where the minus sign comes from commutation of fermionic fields. So we can write 
in the general case 

exp[SA{ip,^)] = l+SA{iJ,ij), (B.5) 
where SaH^^'^P) is defined like SAi'ip^'ip) but with the parameter w'^ replaced by 

f^A + WA;iWA-j - WA-ijWA-ji li A = {z, j} 

Wa = { (B.6) 

\w\ if |A| > 3 

Consider now a Grassmann integral of the form 



Viij.ij) Oij exp 



AGS 



(B.7) 



where t = (ti)igy are parameters, / = (ii,22, . . . ,4) e and J = (ji, j2, • • • ,jk) e 
V'' are ordered /c-tuples of vertices, and 

Oi,j := ^iiV'ii ■■■^ifeV'ife (B.8) 

[cf. (16.11) ]. Here the ii, . . . ,ik must be all distinct, as must the ji, . . . , jt, but there can 
be overlaps between the sets I = {ii, ^2, • • • , ik} and J = {ji, j2, • • • , jfcjo We intend to 
show that flB.7l) can be interpreted combinatorially as a generating function for rooted 



Please note the distinction between the ordered /c-tuple I — (ji, i2, . . . , ik), here written in itahc 
font, and the unordered set I — {ii, i2, . . . , ife}, here written in sans-serif font. 
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orienteco spanning sub(hyper)graphs of G, in which each connected component is 
either a (hyper-)tree or a (hyper-)unicychc. In the case of a unicychc component, the 
rest of the component is oriented towards the cycle, and no vertex from lUJ hes in the 
component. In the case of a tree component, either (a) no vertex from I U J is in the 
component, and then there is either a special "root" vertex or a "root" hyperedge, all 
the rest of the tree being oriented towards it, or (b) the component contains a single 
vertex from I fl J, which is the root vertex, and the tree is again oriented towards it, 
or (c) the component contains exactly one vertex from I and one from J, a special 
oriented path connecting them, and all the rest is oriented towards the path. The 
weight of each configuration is essentially the product of tj for each root i ^ I U J and 
an appropriate weight {w% WA;i or wa-ij) for each occupied hyperedge, along with a 
— sign for each unicyclic using the WA-i/s and a single extra ± sign corresponding to 
the pairing of vertices of I to vertices of J induced by being in the same component. 
(This same sign appeared already in Section [6l) 

Kirchhoff 's matrix-tree theorem arises when all the hyperedges A have cardinality 
2 (i.e. G is an ordinary graph), I = J = {i^} for some vertex iQ, all ti = 0, all = 0, 
and WA-i = WA;ij = wa- The principal-minors matrix-tree theorem is obtained by 
allowing I = J of arbitrary cardinality k, while the all-minors matrix-tree theorem is 
obtained by allowing also I 7^ J. Rooted forests with root weights ti can be obtained by 
allowing ti 7^ 0. On the other hand, unrooted forests are obtained by taking all ti = A, 
I = J = 0, w\ = —Xwa and the rest as above. [More generally, unrooted hyperforests 
are obtained by taking all tj = A, I = J = 0, w*^ = —X{\A\ — 1)wa and the rest as 
above.] The sequences / and J are used mainly in order to obtain expectation values 
of certain connectivity patterns in the relevant ensemble of spanning subgraphs. 

Let us now prove all these statements, and give precise expressions for the weights 
of the configurations, which until now have been left deliberately vague in order not 
to overwhelm the reader. 

We start by manipulating ( IB. 71) . exponentiating the action to obtain 



or, expanding the last products. 



H'^A , (B.IO) 



E'CE 



where I UV consists of the sequence / followed by the list of elements of V in any 
chosen order, and J UV consists of the sequence J followed by the list of elements 
of V in the same order. 



We shall define later what we mean by "orienting" a hyperedge A: it will correspond to selecting 
a single vertex i G A as the "outgoing" vertex. 
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Factors coming from O^jv 



' JUV 



root vertex 
sink vertex 
source vertex 



Factors coming 


from Yl'^A 






root hyperedge 


TA-^i 




pointing hyperedge 


Ipi'^^jTA-^liJ} 




dashed hyperedge 



Table 1: Graphical representation of the various factors in the expansion (IB.lOp . 



We now give a graphical representation and a fancy name to each kind of monomial 
in the expansion (IB.IOI) . as shown in Table [H Please note that in this graphical rep- 
resentation a solid circle • corresponds to a factor ipiipi, an open circle o corresponds 
to a factor ipi, and a cross x corresponds to a factor tpi. 

According to the rules of Grassmann algebra and Grassmann-Berezin integration, 
we must have in total exactly one factor ipi and one factor ipi for each vertex i. 
Graphically this means that at each vertex we must have either a single • or else the 
superposed pair ® (please note that in many drawings we actually draw the o and x 
slightly split, in order to highlight which variable comes from which factor). At each 
vertex i we can have an arbitrary number of "pointing hyperedges" pointing towards 
i, as they do not carry any fermionic field: 




Aside from pointing hyperedges, we must be, at each vertex i, in one of the following 
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situations (Figure [2]): 



1. If 2 G V or i G In J [resp. cases (a) and (b) in the figure], the quantity C/uy',JuF' 
provides aheady a factor ipi4>i'i therefore, no other factors of ipi or ipi should come 



2. If 2 G I \ J, the quantity C/uv'.Juv provides already a factor ifji] therefore, the 
expansion of Y\Sa must provide ipi^ i.e. we must have one dashed hyperedge 
pointing from i. 

3. If z G J \ I, the quantity O/uv.Juv provides already a factor ipi, therefore, the 
expansion of Y\ Sa must provide ■j/'j, i.e. we must have one dashed hyperedge 
pointing towards i. 

4. If z ^ lUJUl^', then the quantity C/uv'.Juv" provides neither ifji nor ipi, therefore, 
the expansion of Y\ Sa must provide both ipi and ipi, so that at i we must have 
one of the following configurations: 

a) a non-pointed vertex of a pointing hyperedge; 

b) a vertex of a dashed hyperedge that is neither of the two endpoints of the 
dashed arrow; 

c) a vertex of a root hyperedge; 

d) two dashed hyperedges, one with the arrow incoming, one outgoing. 

Having given the local description of the possible configurations at each vertex z, 
let us now describe the possible global configurations. Note first that at each vertex 
we can have at most two incident dashed arrows, and if there are two such arrows then 
they must have opposite orientations. As a consequence, we see that dashed arrows 
must either form cycles, or else form open paths connecting a source vertex of I \ J 
to a sink vertex of J \ I. Let us use the term root structures to denote root vertices, 
root hyperedges, cycles of dashed hyperedges, and open paths of dashed hyperedges. 

As for the solid arrows in the pointing hyperedges, the reasoning is as follows: 
If a pointing hyperedge A points towards i, then either i is part of a root structure 
as described above, or else it is a non-pointed vertex of another pointing hyperedge 
'■p{A). We can follow this map iteratively, i.e. go to (p{ip{A)), and so on: 



Because of the finiteness of the graph, either we ultimately reach a root structure, 
or we enter a cycle. Cycles of the "dynamics" induced by correspond to cycles of 



from the expansion of H'^A- 
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Figure 2: Possible ways of saturating the Grassmann fields on vertex i (indicated by 
the small gray disk). 



the pointing hyperedges. We now also include such cycles of pointing hyperedges as 
a fifth type of root structure (see Figure [3] for the complete list of root structures). 

All the rest is composed of pointing hj^eredges, which form directed arbores- 
cences, rooted on the vertices of the root structures. In conclusion, therefore, the 
most general configuration consists of a bunch of disjoint root structures, and a set 
of directed arborescences (possibly reduced to a single vertex) rooted at its vertices, 
such that the whole is a spanning subhypergraph H of G. 

As each root structure is either a single vertex, a single hyperedge, a (hyper-)path 
or a (hyper-)cycle, we see that each connected component of H is either a (hyper-)tree 
or a (hyper-)unicyclic. Furthermore, all vertices in I U J are in the tree components, 
and each tree contains either one vertex from I and one from J (possibly coincident) 
or else no vertices at all from I U J. 

We still need to understand the weights associated to the allowed configurations. 
Clearly, we have a factor WA-i per pointing hyperedge in the arborescence. Root 
vertices coming from V have factors tj, and root hyperedges have factors w\. Cycles 
7 = {io, Ai,ii, A2, = io) of the dynamics of ip {bosonic cycles) have a weight 

u!Ai;ii ■ ■ ■WAi;ii- AH the foregoing objects contain Grassmann variables only in the 
combination and hence are commutative. Finally, we must consider the dashed 
hyperedges, which contain "unpaired fermions" ipi and ipj, and hence will give rise 
to signs coming from anticommutativity. Let us first consider the dashed cycles 
7 = (io, ^1, "^1, ^2? ••••5 H = io)) and note what happens when reordering the fermionic 
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® 

root vertex 


root hyperedge 


source-to-sink path of dashed hyperedges 


• 

cycle of dashed hyperedges 


A = (f)^{A) 
cycle of pointing hyperedges 



Figure 3: The five kinds of root structures. 



fields: 

because ■^j^ had to pass through 21—1 fermionic fields to reach its final location. This 
is pretty much the result one would have expected, hut wc have an overall minus sign, 
irrespective of the length of the cycle (or its parity), which is in a sense "non-local", 
due to the fermionic nature of the fields and ip. For this reason we call a dashed 
cycle a fermionic cycle. 

A similar mechanism arises for the open paths of dashed hyperedges 7 = (zq, Ai, ii, 
A2, . . . where io is the source vertex and i£ is the sink vertex. Here the weight 
WAi;iQhWAr,hi2 ■ ■ ■ multiplies the monomial iJ^,i'i-,iJi-,i'i.,iJi., ■ ■ ■ i'i,_,iJi,_,i>i,, in 

which the only unpaired fcrmions are ipi^ and ipi^. in this order. Now the monomials 
for the open paths must be multiplied by and each source (resp. sink) vertex 
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from an open path must correspond to a vertex of I (resp. J). This pairing thus 
induces a permutation of {1, . . . , k}, where k = |l| = |J|: namely, v is connected by 
an open path to j,r(r) • We then have 

' k \ / ^ \ 

^r=l / \r=l / 

where the first product is Ojj and the second product comes from the open paths. 
This can easily be rewritten as 

k k 

n^V^i.V'v^i.M = Yli'iri^iri'Mr)'Pjr (B.13a) 

-1 

tl^^r^^^ ifli^M.^^j] (B.13b) 



r=l r=l 

k 



\r=l I \r=l 



\r=\ I \r=l 



Putting everything together, we see that the Grassmann integral (ETj) can be 
represented as a sum over rooted oriented spanning subhypergraphs H of G, as follows: 

• Each connected component of H (the unoriented subhypergraph corresponding 
to is either a (hyper-)tree or a (hyper-)unicyclic. 

• Each (hyper-)tree component contains either one vertex from I (the source ver- 
tex) and one from J (the sink vertex, which is allowed to coincide with the 
source vertex), or else no vertex from I U J. In the latter case, we choose either 
one vertex of the component to be the root vertex, or else one hyperedge of the 
component to be the root hyperedge. 

• Each unicyclic component contains no vertex from I U J. As a unicyclic, it 
necessarily has the form of a single (hyper-)cycle together with (hyper-)trees 
(possibly reduced to a single vertex) rooted at the vertices of the (hyper-)cycle. 

• Each hyperedge other than a root hyperedge is oriented by designating a vertex 
i{A) E A as the outgoing vertex. These orientations must satisfy following rules: 

(i) each (hyper-)tree component is directed towards the sink vertex, root ver- 
tex or root hyperedge, 

(ii) each (hyper-)tree belonging to a unicyclic component is oriented towards 
the cycle, and 

(iii) the (hyper-) cycle of each unicyclic component is oriented consistently. 
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Thus, in each (hyper-)tree component the orientations are fixed uniquely, while 
in each unicyclic component we sum over the two consistent orientations of the 
cycle. 

The weight of a configuration H is the product of the weights of its connected com- 
ponents, which are in turn defined as the product of the following factors: 

• Each root vertex i gets a factor ti. 

• Each root hyperedge A gets a factor 

• Each hyperedge A belonging to the (unique) path from a source vertex to a 
sink vertex gets a factor WA;ij, where j is the outgoing vertex of A and i is the 
outgoing vertex of the preceding hyperedge along the path (or the source vertex 
if A is the first hyperedge of the path) . 

• Each hyperedge A that does not belong to a source-sink path or to a cycle gets 
a factor WA;i{A) [recall that i{A) is the outgoing vertex of ^4]. 

• Each oriented cycle (io, Ai, ii, A2, it = io) gets a weight 

e £ 
Yl^A^-i^ - Y[wA^;i„_^i„ . (B.14) 

a=l a=l 

• There is an overall factor sgn(7r). 

B.2 Special cases 

The contribution from unicyclic components cancels out whenever nl=i '^Aa;ia — 
na=i for every oriented cycle (^0,^1,^1,^2, ....,it = io). In particular, this 
happens if WA;ij = wa-j for all A and all i,j G A. More generally, it happens if wa-ij = 
WA;j G^p{<pA;ij) where has "zero circulation" in the sense that Yla=i4>Aa;ia-iia — 
for every oriented cycle (io, Ai, ii, A2, = io). Physically, (p can be thought of as 
a kind of "gauge field' to which the fermions are coupled; the zero-circulation 
condition means that is gauge-equivalent to zero. Note, finally, that if wa-iWa-j — 
WA;ijWA:ji for all i,j e A, then w\ = w\. 

At the other extreme, if we take all tj = 0, all w\ = and I = J = 0, then all 
tree components disappear, and we are left with only unicyclics. 

In certain "symmetric" circumstances, we can combine the contributions from tree 
components having the same set of (unoriented) hyperedges but different roots, and 
obtain reasonably simple expressions. In particular, suppose that the weights wa-i 
are independent of i (let us call them simply wa) , and consider a tree component T 
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that does not contain any vertices of I U J. Then we can sum over all choices of root 
vertex or root hyperedge, and obtain the weight 



A further simplification occurs in two cases: 

• If all ti = t and all = 0, then the second factor in flB.lSp becomes simply 
i|l^(T)|: we obtain forests of vertex-weighted trees. 

• If all ti = t and w\ = t{l — \A\)wa for all A, then the second factor in flB.lSp 
becomes simply t (by virtue of Proposition l2.1l) and we obtain unrooted forests. 

Recall, finally, that if we also take WA-ij = wa for all A and all i,j G A, then the 
unicyclic components cancel and w*a = Wa, so that (1B.15I) reduces to (15.91) . 

It is instructive to consider the special case in which G is an ordinary graph, 
i.e. each hyperedge A G -E is of cardinality 2. If we further take all w\ = 0, then 
the quantity in the exponential of the functional integral flB.7p is a quadratic form 
Slipjip) + = ipMip, with matrix 



Our result for / = J = then corresponds to the "two-matrix matrix-tree theorem" 
of Moon [5, Theorem 2.1] with rik = W{i^k};k for i 7^ k, ra = 0, Sij = W{ijj-ij for i j 
and Sii = 
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